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SECTION  I 


INTRODUCTION 

In  1976  Montrose,  Moynihan  and  Sasabe1  presented  a  model 
for  the  behavior  of  a  viscoelastic  fluid  entrained  in  a  concen¬ 
trated  line  contact.  Their  calculation  gave  a  result  for  the 
traction  coefficient  (the  ratio  of  the  tractive  to  the  normal 
force)  at  low  slip  in  terms  of  an  average  effective  viscosity  {T. 
Assuming  a  Hertzian  semi-elliptic  pressure  profile  in  the  EHD 
contact  their  result  can  be  written  in  the  form 

CT  =  (8R/nhE)  [2i7/t0  J  (AU/U).  (1) 

Here  CT  is  the  traction  coefficient,  R  =  Rj  R2  /  ( Rj  +  R2 )  is  an 
effective  radius  for  the  disk  pair  (simple  twin  disk  geometry 
was  assumed)  ,  of  radii  RJ  and  R2  ,  E  is  an  effective  elastic 
modulus  for  the  disks  (E  is  related  to  Young's  modulus  EQ  by  E  = 
Eq(1  -  v2)  where  v  is  Poisson's  ratio),  tQ  is  the  transit  time 
of  a  fluid  element  through  the  contact,  the  thickness  of  which 
is  2h.  The  average  velocity  of  the  disks  is  U,  i.e.,  U  =  (U;  + 
U2 )/2  where  Uj  and  U2  are  the  individual  disk  velocities,  and  AU 
*  Uj  -  U2  is  the  slip.  The  thrust  of  the  Montrose,  Moynihan  and 
Sasabe  paper  was  to  determine  what  properties  of  a  lubricant  are 
instrumental  in  determining  the  average  effective  viscosity  u 
and  thus  the  traction  coefficient.  With  this  kind  of 
information  one  could  (at  least  in  principle)  molecularly 
engineer  a  lubricant  with  viscoelastic  characteristics  that 
would  optimize  f  and  concommitantly  CT.  The  details  of  the 


calculation  do  not  concern  us  here  (indeed  a  recent  revision  of 
the  model  by  Heyes  and  Montrose  in  which  the  result  is 
expressed  entirely  in  terms  of  directly  measureable  quantities 
will  be  used  in  this  paper  for  the  calculations);  rather,  we 
shall  make  use  of  only  a  few  of  their  general  conclusions. 

As  the  lubricant  enters  the  contact  zone  it  is  subjected 
suddenly  to  a  large  pressure  increase  causing  structural  changes 
in  the  liquid;  these  lead  to  a  large  increase  in  the  charac¬ 
teristic  liquid  relaxation  times  of  the  fluid  and  thus,  on  the 
time  scale  of  the  process  (transit  times  typically  ~  10"4  s) , 
the  lubricant  acts  like  a  solid  shear  stress  transmitting  pad. 
If  one  could  presume  that  this  liquid-to-glass  transition 
occurred  instantaneously  upon  the  fluid's  entering  the  contact 
the  average  effective  viscosity  would  be  given  by 

u  -  G«t0/2,  (2) 

where  G,,,  is  the  instantaneous  shear  rigidity  of  the  lubricant. 
From  the  standpoint  of  maximizing  traction  this  is  the  most 
desirable  situation  and  immediately  leads  to  two  general  lubri¬ 
cant  design  principles: 

(1)  The  shear  rigidity  of  the  fluid  G*  should  be 
as  large  as  possible.  Since  G,,,  increases 
with  pressure  as  the  liquid  traverses  the 
contact,  a  large  value  of  the  pressure  co¬ 
efficient  3 G—/3  p  is  also  desirable. 

(2)  The  fluid  should  undergo  the  liquid-to-glass 
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transition  as  rapidly  as  possible  after 
entering  the  contact  zone.  There  are  three 
viscoelastic  parameters  that  primarily 
control  the  rate  at  which  this  transition 
occurs:  these  are  (a)  the  pressure-viscosity 
coefficient,  a  =  3 (In  n )/3 P;  (b)  the  ratio  of 
instantaneous  to  equilibrium  compressibili¬ 
ties  <»A0  ;  and  (c)  the  form  of  the 
viscoelastic  relaxation  function  $(t). 

There  are,  of  course,  several  other  factors  that  influence 
the  behavior.  From  the  viewpoint  of  simply  maximizing  iT,  one 
would  want  the  inlet  shear  viscosity  n0  to  be  large;  however, 
insofar  as  nQ  is  of  major  importance  in  determining  the  film 
thickness,  one  is  limited  to  a  rather  small  range  of  feasible  n0 
values.  The  ratio  of  structural  to  shear  relaxation  times  is 
also  one  of  the  factors  that  determines  iT.  However,  for  the 
type  of  organic  fluids  that  satisfy  the  other  criteria  of 
acceptable  traction  performance,  one  generally  finds  this  ratio 
to  lie  in  the  range  of  about  one  to  ten.  Attempting  to  design 
or  find  a  material  for  which  the  ratio  is  less  than  10“2  or  so 
(which  would  increase  traction  somewhat)  is  just  not  a  viable 
alternative. 

Consequently  our  aim  in  this  research  has  been  to 
characterize  a  set  of  liquids  with  respect  to  the  four 
properties  ( G,,, ,  a,  4>  ,  and  described  above. 
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SECTION  II 


THE  VISCOELASTIC  PARAMETERS  AND  THEIR  MEASUREMENT 

A.  Shear  Rigidity,  Gm . 

For  a  given  set  of  conditions  (rolling  speed,  slip-roll 
ratio,  Hertzian  pressure,  etc. )  the  traction  coefficient  is 
essentially  proportional  to  G,,,  ,  and  consequently  it  is  desirable 
to  have  as  large  a  value  cf  G,,  (and  3Go»/8P)  as  possible. 

To  measure  G,,,  we  take  advantage  of  the  expression 

G.  =  pVTf0„2,  (3) 

giving  the  shear  modulus  in  terms  of  the  density  p  and  the 
velocity  of  propagation  of  a  transverse  wave  vt,<»  '  under 
conditions  for  which  the  period  of  the  wave  is  much  shorter  than 
the  shear  relaxation  time. 

High  frequency  wave  speeds  (both  transverse  and  longitu¬ 
dinal)  were  measured  using  the  technique  of  Brillouin 
scattering.  Here  one  illuminates  the  sample  with  a  collimated 
beam  of  monochromatic  light  and  observes  the  spectrum  of  the 
light  scattered  through  some  angle  0  .  This  spectrum  consists  of 
a  triplet  symmetrically  placed  about  the  incident  light  fre¬ 
quency,  u>0  ;  the  two  shifted  Brillouin  lines  appear  at  fre¬ 
quencies  u) Q  ±  u>L  if  the  incident  polarization  is  normal  to  the 
scattering  plane.  The  Brillouin  frequencies  <oT  and  give, 
respectively,  the  frequencies  of  transverse  and  longitudinal 
thermal  phonons  with  a  wavevector  k  satisfying  the  condition 
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k 


(  4n n/X )  sin(0 /2 ) 


(4) 


where  n  is  the  refractive  index  of  the  material  under  investi¬ 
gation  and  X  is  the  wavelength  of  the  incident  light.  Thus  from 
the  measured  Brillouin  shifts  one  obtains  the  transverse  wave 
speed 

Vfp  =  gj  T/k  =  ( G)  fpX  /4u  n)/sin(0  /2  )  (5) 

and  the  longitudinal  wave  speed 

VL  =  wj/k  =  (wLX/4irn)/sin(0/2)  .  (6) 

At  sufficiently  high  pressures  and/or  low  temperatures  where  the 
viscoelastic  relaxation  times  are  long  (>  10“9  s)  ,  these 

approach  the  limiting  wave  speeds  VT f<a  and  VLf00. 

In  our  measurements  an  argon  ion  laser  operating  at  X  = 
514.5  nm  was  used  as  the  source  and  a  piezoelectrically  scanned 
Fabry-Perot  interferometer  was  used  as  the  spectrum  analyzer. 
The  scattering  angle  was  90°.  Refractive  indices  were  measured 
using  a  modified  Abbe-Grauer  ref ractometer .  These  measurements 
were  extrapolated  to  high  pressures  using  the  Lorenz-Lorentz 
formula 

(n2  -  1)  /  (n2  +  2)  =  YP  , 

I 

.  where  y  is  a  constant.  The  extrapolated  indices  agreed  to 

within  about  0.1%  with  values  obtained  by  measuring  the 
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deflection  of  the  laser  beam  when  a  small  prism  was  in  the  high 
pressure  cell.  The  values  of  V.^,.  obtained  are  estimated  to  be 
accurate  to  about  ±  1%.  Densities  were  measured  at  atmospheric 
pressure  to  an  accuracy  of  ±  0.1%  using  a  pycnometer.  At 
elevated  pressures  a  densitometer  essentially  similar  to  that 
described  by  McDuffie  et  al.^  was  used  to  obtain  values  accurate 
to  ±  0.2%.  The  overall  precision  of  the  measured-shear  moduli 
is  better  than  ±  3%. 

B.  Pressure-Viscosity  Coefficient,  a 

Because  large  traction  coefficients  are  obtained  when  the 
lubricant  becomes  glassy  soon  after  entering  the  contact  zone, 
the  pressure  variation  of  the  viscosity  is  crucial  in 
determining  the  lubricant's  traction  characteristics.  For  the 
materials  being  reported  here  (and  for  many  other  lubricant 
fluid  candidates)  ,  the  pressure  dependence  of  the  viscosity  is 
described  reasonably  well  by  the  simple  Barus  formula 

n  =  A  exp(aP) , 

where  A  is  a  constant  and  a  =  3 (In  n )/3P  is  called  the  pressure- 
viscosity  coefficient.  Shear  viscosities  in  this  work  were 
measured  at  atmospheric  pressure  to  an  accuracy  of  ±  0.1  in 
log  n  using  a  Cannon-Fenske  viscometer.  At  higher  pressures  a 
"falling  slug"  viscometer  described  by  McDuffie  and  Barr-*  was 
employed . 
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C .  Viscoelastic  Relaxation  Function,  <t>  (t) 

The  dynamic  response  of  the  lubricant  to  shearing  and 
compression  is  described  by  time-dependent  relaxation  functions. 
In  the  simplest  case  the  relaxation  functions  take  the  form  of  a 
simple  exponential  decay: 


<|>(t)  =  exp(-t/i  ) . 


Usually  this  simple  form  is  inadequate  to  describe  the  behavior 
and  a  sum  or  continuous  distribution  of  such  exponentials  is 
used : 


(t)  =  Ig i  expf-t/ii)  -*•  /“dx  g(x  )  exp(-t/x  ) . 
i  0 


The  width  of  the  distribution  function  g(x)  is  an  important 
factor  in  determining  the  traction  coefficient — the  narrower  is 
g(x),  i_.e.  ,  the  closer  <t>(t)  is  to  an  exponential  decay,  the 
larger  will  be  CT.  A  convenient  parameter  characterizing  this 
width  is  the  variance  B  defined  by 


B  =  [  <x2  >  -  <x>2]1/2<x>, 


where 


/  dt  <Mt) 
o 


/  dx  g(x  )  x  ; 
o 


(10a) 


<xz  > 


=  /"dt  <{>  (t)  =  /"dx  g(x  )  x2 
0  0 


(10b) 
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For  most  traction  lubricants  B 

ranges  from 

about 

one 

to  three; 

values  substantially  less 

than  unity  will 

lead 

to 

unusually 

large  traction  coefficients. 

For  the  materials 

studied 

in  this 

work 

the 

measured 

relaxation  functions  were 

found 

to  fit  rather 

well 

to  the 

empirical  form- 


4>(t)  =  exp[ -( t/x  )B] , 


(ID 


where  0  <  0  <  1.  In  terms  of  the  parameters  tq  and  6  one  has 
<t>  =  (t0/b)  r  (1/0  )  ,  (12a) 

and 

<t2>  =  (t2/0 )  r (1/0  ).  (12b) 

The  relaxation  width  B  is  given  by 

b  =  [ 0 r ( 2/0  )  -  r (I/0 )2 ] 1/2  /  r (I/0 ) .  (13) 

In  Eqs.  (12)  and  (13),  r(...)  denotes  the  gamma  function. 

The  measurement  of  the  relaxation  function  $(t)  was  carried 
out  using  the  technique  of  digital  correlation  spectroscopy.6-9 
Here  using  a  digital  correlator  one  measures  a  function  F(t) 
describing  the  correlation  of  intensity  fluctuations  of  light 
scattered  by  the  liquid  under  investigation.  The  directly 
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measured  function  F(t)  is  related  to  the  viscoelastic  relaxation 
function  <{>{t)  by  the  relationship 

F(t)  =  a  +  b  <j>(t)2,  (14) 

where  a  and  b  are  constants  determined  by  the  experimental 
conditions . 

D.  Instantaneous-Equilibrium  Compressibility  Ratio  k^A  . 
The  ratio  ic^Ag  determines  the  extent  to  which  the  lubri¬ 
cant's  structure  changes  "instantaneously"  in  response  to  the 
pressure  impulse  encountered  as  it  passes  through  the  EHD 
contact  zone.  The  largest  possible  traction  coefficient  would 
result  if  the  entire  response  were  instantaneous  allowing  the 
viscosity  to  increase  more  rapidly;  this  behavior  would  occur  in 
the  limit  *„A0  ♦  1.  Typically  for  fluid  lubricants  one  finds  a 
value  of  about  0.5  to  0.6  to  be  more  common. 

The  equilibrium  compressibility  was  obtained  by  carrying 
out  ultrasonic  "sing-around"  measurements  of  the  low  frequency 
sound  speed  and  then  employing  the  relationship 

<0  =  1  /  PVQ2.  (15) 

The  instantaneous  compressibility  was  obtained  from  the 

formula 

K.  =  [M»  -  4G„/3]-1  ,  (16) 

where  Mw ,  the  limiting  longitudinal  modulus,  was  obtained  from 

/  ^ 


9 


the  longitudinal  wave  speed  VL>U>  measured  using  Brillouin 
scattering  ( =  pVLf-2).  The  errors  in  the  determination  of  kq 
(±  2%),  M,,,  (±  2%)  and  G*  (±  3%)  combine  so  that  the  ratio  k„A0 
can  be  judged  to  be  accurate  to  about  ±  6%. 
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SECTION  III 


EXPERIMENTAL  RESULTS,  ANALYSIS  AND  INTERPRETATION 

A.  Materials 

Three  lubricant  materials  were  evaluated  for  their  traction 
characteristics.  They  were  designated  as  (a)  MLO-78-187,  a 
MIL- L-7808  base  stock,  (b)  MLO-77-127,  a  polyalphaolef in,  and 
(c)  MLO-76-136a,  a  narrow  molecular  weight  (574  avg)  petroleum 
fraction.  These  materials,  provided  by  the  Air  Force  Materials 
Laboratory,  Fluids,  Lubricants  and  Elastomers  Branch,  were 
filtered  through  a  Millipore  filter  (pore  size  0.65  uni)  to 
remove  dust  particles  and  were  then  used  without  further 
purification. 

B.  Presentation  of  the  Results 

Longitudinal  modulus  data  were  obtained  at  atmospheric 
pressure  as  a  function  of  temperature  and  at  room  temperature  as 
a  function  of  pressure  for  each  of  the  three  sample  materials. 
These  data  are  presented  in  Figs.  1  and  2. 

Shear  modulus  data  could  be  obtained  only  for  two  of  the 
three  liquids  (MLO-77-127  and  MLO-76-136a)  and  only  at  the 
highest  pressures  studied.  Lower  pressure  values  were 
extrapolated  by  assuming  that  the  ratio  —  MQ  )  is 
independent  of  pressure  as  is  the  case  for  many  liquids.  For 
the  third  liquid,  MLO-78-187,  the  shear  modulus  was  estimated  by 
assuming  that  the  ratio  -  MQ  )  =  0.35  as  was  found  for 
MLO-77-127.  This  value  is  typical  of  that  found  in  many 
liquids.  The  measured,  extrapolated  and  estimated  shear  modulus 
data  are  given  in  Table  I. 
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Density  measurements  versus  temperature  and  pressure  for 
the  three  liquids  are  shown  in  Fig.  3.  These  data  were  used  in 
the  computation  of  the  moduli  and  in  the  extrapolation  of  the 
refractive  indices  according  to  the  Lorenz-Lorentz  formula.  The 
relevant  data  for  this  extrapolation  are  given  in  Table  2. 

Shear  viscosity  data  versus  temperature  and  pressure  are 
shown  in  Figs.  4  and  5.  It  is  evident  from  the  latter  figure 
that  over  the  pressure  range  of  the  measurements  the  data  can  be 
described  reasonably  well  by  the  Barus  formula. 

The  measurements  of  the  viscoelastic  relaxation  function 
4>(t)  were  carried  out  for  MLO-78-187  and  MLO-77-127  at  low 
temperatures  and  elevated  pressures.  These  data  are  shown  in 
Figs.  6  and  7.  It  was  not  possible  to  obtain  an  accurate 
determination  of  <f>(t)  for  MLO-76-136a  since  the  scattered 
intensity  increased  by  nearly  a  factor  of  ten  over  the  course  of 
an  eight-hour  run.  A  determination  of  the  cause  of  this 
phenomenon  is  outside  the  scope  of  this  program.  An  attempt  was 
made  to  obtain  an  approximate  form  for  <f>(t)  by  carrying  out  a 
"quick"  run  at  21°  C  and  0.48  GPa.  These  data  were  fit  to  the 
fractional  exponential  form  in  Eq.  (11)  with  parameters  tq  =  8  x 
10”5  s  and  p  =  0.4.  In  the  analysis  given  below  the  forms  found 
for  4>  ( t )  were  assumed  to  be  pressure  and  temperature  independent 
and  the  pressure  variation  of  relaxation  times  was  assumed  to  be 
similar  to  the  shear  viscosity  as  is  usually  the  case. 
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C.  Analysis  and  Interpretation 

As  we  have  emphasized  above,  the  objective  of  this  study 
was  to  evaluate,  in  terms  of  fundamental  viscoelastic  data,  the 
traction  characteristics  of  the  candidate  fluids  for  application 
to  rolling  bearing  analysis.  To  establish  a  basis  for  this 
evaluation,  it  is  useful  to  recall  the  results  of 
characterization  studies  similar  to  this  one  that  have  been 
carried  out  on  lubricants  with  large  traction  coefficients.  In 
Table  III  we  summarize  the  appropriate  viscoelastic  data  for 
Santotrac  40  (a  commercial  traction  fluid  manufactured  by  the 
Monsanto  Corporation)  and  the  synthetic  lubricant  5P4E  (five 
phenyl  rings  joined  by  ether  bridges).  We  have  also  included  in 
this  table  the  corresponding  data  for  a  short  chain  methyl- 
phenylsiloxane  polymer  MPS-710.  (Dow  Corning  silicone  fluid 
No.  710).  The  slopes  of  the  plots  of  traction  coefficient  CT 
versus  slide-roll  ratio  AU/U  were  computed  in  the  low  slip 
region  using  the  model  of  Heyes  and  Montrose.  For  these 
computations  the  following  set  of  typical  external  parameters 
were  used: 

R  =  1 . 9 1  cm 
E  =  233  GPa 
Pjj2  =  1.0  GPa 

The  fact  that  the  traction  slope  of  5P4E  is  about  30%  higher 
than  that  of  Santotrac  40  is  in  rough  agreement  with  what  is 
experimentally  observed.  The  larger  value  obtained  for  5P4E  is 
primarily  attributable  to  the  smaller  value  of  B  and  the  larger 
value  of  <(>,/k0.  The  smaller  value  of  G.  for  5P4E  causes  a 
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reduction  in  the  difference  of  the  two  traction  slopes  (if  equal 
values  of  G„  are  assumed,  the  traction  is  50%  larger  for  5P4E 
than  for  Santotrac  40). 

For  the  three  materials  characterized  in  this  report,  the 
traction  slopes  are  less  than  half  of  those  found  in  5P4E  and 
Santotrac  40.  The  origins  of  these  lower  values  are  rather 
easily  understood  in  terms  of  a  comparison  with  Santotrac  40. 
Consider  first  MLO-78-187.  Its  shear  modulus  is  1.1  GPa 
compared  with  1.9  GPa  for  Santotrac  40.  This  alone  would  effect 
a  reduction  of  the  computed  traction  slope  from  102  to  61.  The 
very  small  pressure  viscosity  coefficient  for  MLO-78-187  (15 
GPa-1  vs.  40  GPa-1  for  Santotrac  40)  leads  also  to  a  roughly 
50%  reduction  in  the  traction  slope  as  does  the  very  broad 
relaxation  function  width.  For  MLO-77-127  and  MLO-76-136a  the 
shear  moduli  are  somewhat,  although  not  drastically,  lower  than 
that  for  Santotrac  40  and  this  leads  to  some  reduction  in  the 
traction. 
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SECTION  IV 


SUMMARY 

The  viscoelastic  properties  of  the  three  potential  traction 
fluid  candidates  studied — MLO-78-187,  MLO-77-127,  and 
MLO-76-136a  — have  been  found  to  predict  traction  coefficients 
somewhat  less  than  those  of  the  best  traction  lubricants  now 
available.  The  principal  cause  of  the  lower  traction  coeffi¬ 
cients  is  the  rather  small  pressure  dependence  of  the  viscosity 
for  these  fluids;  in  addition,  their  measured  shear  rigidity 
moduli  are  slightly  less  than  those  in  fluids  specifically 
designed  to  achieve  high  traction. 


SECTION  V 


TABLES 


TABLE  1 


SHEAR  MODULI  FOR  THE  THREE  CANDIDATE  LIQUIDS  UNDER  INVESTIGATION 


Pressure 

MLO-78-187 

Shear  Modulus  (GPa) 
MLO-77-127 

MLO-76-136a 

( GPa) 

(21°  C) 

(22°  C) 

(21°  C) 

o 

• 

tvj 

0.7* 

1 . 4  + 

1. 7+ 

0.4 

1.1* 

2.0 

2.2 

♦Estimated  by  taking  -  MQ  )  =  0.35 


♦Extrapolated  by  assuming  that  G./fM*  -  Mn )  is  independent 
of  pressure.  u 
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TABLE  2 


REFRACTIVE  INDEX  DATA  FOR  THE  THREE  CANDIDATE  LIQUIDS  STUDIED 

MLO-78-187  MLO-77-127  MLO-76-136a 


Refractive  index 

(20°  C,  1  atm)  1.48661  1.49301  1.48975 

Lorenz-Lorentz 
coefficient,  y 

(cm3/g )  0.3090  0.3491  0.3371 


TABLE  3 


VISCOELASTIC 

PARAMETERS  AND 

COMPUTED 

TRACTION 

SLOPES 

FOR 

CANDIDATE 

LUBRICANTS 

(a) 

G. 

(GPa) 

a 

(GPa-1 

B 

) 

*  00  /*  Q 

Traction 

slope 

Santotrac  40(*')) 

1.9 

40 

1.33 

0.5 

102 

5P4E  < c ) 

1.6 

45 

0.88 

0.74 

129 

MPS-710(d) 

1.4 

45 

>2.5 

0.52 

<45 

MLO-78-187 

0.7 

14 

7.5 

0.66 

18 

MLO-77-127 

1.4 

16.5 

1.4 

0.58 

48 

MLO-76-136a 

1.7 

23 

2.1 

0.41 

50 

(a)  Parameters  given  are  for  20°  C,  0.2  GPa.  Where  data  was 
not  specifically  available  at  this  pressure  it  was  obtained 
by  extrapolation  from  higher  or  lower  pressure. 

(b)  Data  from  R.  Moeller,  R.  Meister,  and  C.  J.  Montrose, 

Technical  Report  #8,  ONR  Contract  No.  N0017-75-C-0585 
(  1978)  . 

(c)  Data  from  J.  F.  Dill,  P.  W.  Drake,  and  T.  A.  Litovitz, 

Trans.  ASLE  1_8,  209  (1975);  Tech.  Reports  Nos.  2  and  3, 
ONR  Contract  N00017-67-A-0377-0018 . 

(d)  Data  from  B.  Dorn,  Ph.  D.  thesis.  Catholic  University 

(  1977).  The  form  of  <f>  ( t )  is  found  to  be  a  superposition  of 
that  given  in  Eq.  (11)  with  e  s  0.40  plus  a  long  time  tail 
varying  as  t-rt  with  n  =  0.25.  This  latter  feature  was 
neglected  in  the  computation  of  B  and  of  the  traction 
slope.  The  traction  slope  thus  obtained  represents  an 
upper  bound. 


19 


SECTION  VI 


ILLUSTRATIONS 
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MODULUS,  GPa 


MODULUS,  GPa 


FIGURE  7.  The  measured  relaxation  function  <Mt)  is  shown  for  MLO-77-127  for 
T  =  -16 °C  and  P  =  0.48  GPa.  The  solid  curve  is  computed_from 
Eq.  (11)  using  the  parameters  8=  0.5  and  r  =  9.5  10  l*s. 
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